Singular three-body integrals for Hylleraas wave functions 
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Abstract 
In this paper, we study singular integrals that appear in calculations of relativistic and quantum 
electrodynamic corrections for three-body Coulomb systems using Hylleraas coordinates. A general 


scheme is developed to separate the divergent part of a singular integral involving Hylleraas basis 


functions of type rirjrj. exp(—ar1 — Br2). The integrals of type Ip(a, b, —2 — n; a, 8) with a > —1, 


b > —1, and n = 1, 2,3 are considered and matched with T; m (a, 8,7) for both singular and regular 
parts. Moreover, we also consider the singular integrals of type I,(a, b, c; œ, 8) where c > —1 and 


at least one of a and b is negative. 


* zxzhong@wipm.ac.cn 


I. INTRODUCTION 


Few-body atomic and molecular systems are basic grounds for testing physical laws, such 
as quantum electrodynamics (QED), as well as are one of the sources of determining funda- 
mental constants of nature by combining precision measurements and precision theoretical 
calculations. Important examples include the atomic helium whose 2? P; fine structure holds 
potential for deriving the fine structure constant a [1] 2], the hydrogen molecular ions Hy 
and HDt from which one can determine the electron to nuclear mass ratios [3] [4], and 
the antiprotonic helium pHet that can be use to test the CPT theorem [4] in relativistic 
quantum field theory. 

From theoretical side, the eigenvalue problems for the nonrelativistic Hamiltonians of 
these Coulomb three-body systems can be solved very precisely using variationally con- 
structed wave functions in Hylleraas coordinates, which lays a foundation for further evalu- 
ating relativistic and QED corrections to a required precision. The matrix elements of many 
singular relativistic and QED operators can be divergent individually; however, the diver- 
gent parts of these singular integrals should be cancelled with other counterparts, resulting 
in a finite value. In general, rotational degrees of freedom in a matrix element can first be 
isolated analytically and treated using angular momentum theory, and the remaining part 
of the matrix element becomes an integral that involves only radial coordinates, such as rı, 
r2, and rı2 = |r) — | for a two-electron system [5]. Of course, the basic type of radial 
integrals depends on the variational basis functions chosen. 

One of the most popular variational basis sets used for three-body Coulombic systems 
is [6] 

Jew hipaa Me ram) (1) 
where Vee (AL, f2) is the angular momentum eigenfunction formed by ¢; and ¢2, and a, b, 
and y are three nonlinear variational parameters. This type of basis set has been successfully 
applied to the calculations of relativistic and QED effects for systems like the antiprotonic 
helium [7] and the hydrogen molecular ions [8]. The radial integrals for this type of basis 
functions are of the form 
o0 oO ry tre 
Pala) = J an | dr n 7 dari rpe ee, (2) 


The above integral is invariant under the simultaneous permutation of l, m, n and a, £6, 


y, and its range of integration covers all possible values of rı, r2, and ry with (r1, r2, 
rı2) forming a triangle. This integral can be evaluated recursively [9], and the analytical 
expressions for the singular integrals P_,o(a, 8, y) and P'_p-1 (a, 8,7) with integer p > 0 
have been derived by Harris et al. [10]. Korobov has recently worked out the singular 
integrals of the type T_2-20(a, 6,7). It should be mentioned that the following recursion 


relation will be used frequently in dealing with singular integrals: 


o 
rt ge p, y) = AAG p, 7) $ (3) 


Another widely used basis set is based on the correlated Hylleraas functions 
ien e, (4) 


which has been applied to high-precision calculations for two-electron atomic systems [12] 
[13], the hydrogen molecular ions [14H18], and the antiprotonic helium [19]. After averaging 
the rotational degrees of freedom, the radial integral can be reduced to 
o0 oO ry tre 

I, (a, b,c; a, B) -f rdn f rədrə J. riod. rororsje oP"? P (cos 0), (5) 
where cos@ = (r? + r? — r?,)/(2rir2) and P,(cos@) is the Legendre polynomial. By using 
one of the recurrence relations for P,(cos@), I,(a,b,c;a,) can be expressed in terms of 
various Io(a,b,c;a,) [5]. The analytical expression for J, has been worked out by Yan 
and Drake for the case where a,b,c > —2 anda+b+c+5 > 0. An analyti- 
cal expression for singular integral Io(a,b, —2 — n;a, 8) with n > 1 was obtained in 
together with some examples for analytical cancellation of some singularities, where the 
singular parts are written in terms of some subsidiary divergent integrals. The analytical 
expression for I9(a,b, —2 — n;a,() has been derived by Drake [13] applying the differen- 
tial chain rule on the generation integral Io(—1,—1,—2 — n;a, 8), where the generation 
integral was treated by replacing the lower limit of the integral by € and studying the be- 
havior as e + 0*. This ideal was adopted by Harris et al. who studied the singular 
integrals of type Timn(a,,y) . For evaluating high-order relativistic and QED effects, 
more singular integrals than Ip(a,b, —2 — n; a, 8) are required, such as I,(—2 — n, b, c; a, B), 
I,(—2 —m, —2 —n,c; a, 8), and I,(—2 — m,b, -2 —n; a, 8) with m,n = 0,1, 2,.... 

Comparing Eqs. and (5), one can establish the following relation 


Iola, b,c; œ, B) = Vesa eer (Q, 8,0). (6) 


Most singular integrals I,(a,b,c;a,2) can be evaluated through Dymn(a, 8,0) except the 
integrals of type I,(—2 — m, —2 — n, c; a, 8) with m,n > 0. For the case of m = n = 0, let 
us evaluate [9(—2, —2, c; a, 8) through the generation function P_1_1(a, 6,0) that can be 
obtained from the following formula (see Appendix F of Ref. [10}) 


_ 1| safety . a+b , a+B “a 
P_1-10(a, 8,7) F f ln (5 J 2) + dilog (=) + dilog (=) fe =| (7) 


by taking the limit y — 0, where dilog is the dilogarithm function. It is, however, impossible 
to evaluate Ip(—2, —2,c;a, 6) with c > 0 by differentiating Io(—2, —2, —1; a, 6) with respect 
to y which equals zero. 

The purpose of this paper is to study two types of singular integrals arising from 
I,(a, b,c; a, 8), including Iola, b, —2 — n;a,8) in Sec. [I| and I,(—2 — m, —2 — n, c; a, 8) 
with c > —1 in Sec. The singular parts of these integrals will be given explicitly. The 


verification of our results will be presented in Sec. [IV] together with a conclusion. 


II. SINGULAR INTEGRALS o(a, b, —2 — n; œ, 8) 


The derivation of Jo(a, b, —2 — n; a, B) with integer n > 1 has been performed in Ref. [12]. 


Here we quote the final result: 


Ip(a, b, —2 — n; a, B) =O, + O2 + O3, (8) 
where 
1(a+1)!(6+1)!\(a+b4+3-—-n)! B-a 
== 2,n;a+b+4; > 
O1 n (a Ae fy cits 3)laat2-ngo+2 2Fı(b 2, n;a b 4; B ) ; p ZQ 
1(a+1){b+1)a+b0+3- n)! on aB 
~ on (a +b+ a ener 2ki(a 2, MaA b+ 4; & ) ’ p <a 


b+1 


i GD (@+b+2—7) 
= n Ly OFI i TA 


Xn- (8) , 
j=0 
a+l1 


i aril Gerta 
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j=0 


In the above, X,(x) is defined as 


which has different analytical forms for different n. Equation is valid for a,b > —1 and 
a+b+3 > n > 1. The divergent part of the integral is identified as the uncalculated 
subsidiary integrals X,(x). In order to simplify singularity cancellation procedure so that 
some matrix elements involving singular operators like 1/r3, and 1/ri, can be evaluated 
more easily, it is desirable to find out the explicit form of X,(z). 

In Ref. [13], Drake presented a different derivation for [o(a,b, —2 — n; a, 3), in which the 
integral is evaluated by differentiating [9(—1, —2, —2 — n; œ, 3) according to 


Fa) a+1 Ə b+1 
Io(a,b, —2 = n; a, B) = (-Ż) (-3) Io(—1, —1, —2 — n; a, B), (10) 


lar e7217 Bro 
Ip(- 1,- — n; a, b) = a ara f ars | — mi ` (11) 
Iri2— a ae) 


The above integral can further be simplified to the following form after performing integra- 


where 


tions over rı and ro 


2 ” d 
Io(—1,—1,—2 — n; a, 8) = — r / oe (eee — ef"), (12) 
an — € 


"2 
where the lower limit of the integral has been replaced by an infinitesimally small positive 
quantity e. One then considers the behavior of the integral as e 0, where the regular part 
of the integral does not depend on € and the singular part depends on € in the form of either 
Ine or 1/e. 
In order to evaluate Eq. (12), it would be convenient to introduce an auxiliary function 


L,(x) defined by 


L(x) = f t Pe "dt. (13) 


For positive integer p, Lp(x) can be recast to the form listed below as € — 0 (see Eq. (9) of 
Ref. and Appendix |A) 


fv) -Inee) (#0) 
Lyla) = (-ayr-t| MES ——_|, p>0 14 
= CA ow! + 245-7 -0) oa 
where w(p) is the digamma function defined by w(p) = —yr + 7”. im -1 with yg = 
0.57721... being the Euler’s constant. For p < 0, we have 
Pepel) 
Lolz) = — Tp ’ pot. (15) 


Using L,(x), the integral in Eq. becomes 


2 


Io(—1,—1,—2 — n; œ, ĝ) = T 


eaaa n>0, (16) 


which is the same as Eq. (11) of Ref. [IO] by taking y as zero. Thus, the expressions for 
Io(—1,—1,—2 — n; a, 8) for the first few values of n can be obtained from Eqs. (13)-(16) of 
Ref. by setting y = 0. Therefore, Tola, b, —2 — n; a, 8) can be calculated according to 


nlab, -2- nia, 8) = (-2) _ (-5) To o Eneila)]- (7) 


It should be mentioned, however, that the integral Tola, b, —2— n; a, 8) can also be evaluated 
recursively using the scheme of Korobov [9]. 

Now we have two sets of expressions for the same singular integral Tola, b, —2 — n; a, 8). 
The first set Eq. can be evaluated directly without use of any recurrence relations; 
however, it is complicated to deal with singularities. On the contrary, the second set Eq. 
can be evaluated recursively because of the derivative nature of the expression, and the 
regular part can be obtained by simply dropping singular terms involving lne and 1/e. We 


can unify these two approaches by choosing the following expressions for X,,(«) that matches 


Eq. with Eq. 


L,(2) p<0 
La n=landp>0 
X,(x) = á (18) 
Latam: n=2andp>0 
L,(2) acon =o 4 (1— ô) n=3andp>0 


The equivalence of Eq. and Eq. can thus be verified explicitly by using a symbolic 
software such as Mathematica. 

It is noted that Eq. will not be applicable for a + b +3 < n due to the existence 
of (a+b+3-— n)! that appears in O,. In this case, Jola, b, —2 — n;a, 8) can be evaluated 
according to Eq. (17), for example, 


Io(—1, -1, —4; a, 8) = asp te — (a + B)((3) — Ind} + zam- amna), a AB 


+ {+—alw(3) —Ine]} +Inat+ 3, a= B 


1 fa e y Seti Ind} + satan B lng — ao’ lna a B 
Io(—1, —1, —5; a, 8) = a+ | e € 3 (4) | sara ( ) À 7 
+ {3 — +a [p(4) —Ine]} — a(3lna +1), a= 8 
(20) 
Other integrals Ip(a, b, —5; a, 3) with a > —1 and b > —1 can easily be evaluated by differ- 
entiating [9(—1, —1, —5; a, 3) with respect to a and 8 repeatedly. 


III. SINGULAR INTEGRALS J,(a,b,c;a,8) WITH c>—1 


The integral I,(a,b,c;a,() may be calculated directly using the method developed by 
Yan and Drake in Ref. [20], in which J,(a, b,c; a, 8) can be reduced to the following form 


Lo 
2 
b, Ceqktr(a, b, c; a, 8; q, k), 21 
nlabas B) = samy Cd eh 0 50,8 e1 
where Cegk and La are respectively 
2+1 f c+2 "E 2k +2- c (22) 
eqk — c+2 2k +1 16 2k +2q—2t+1 
and 
1 g 
=C—q c is even 
Lo =<? (23) 
s(c+1) cisodd, 
and Ip is defined as 
Igla, b, c: a, 3; q, k Di a f dry ef oe rip F q 2ko arı pro (24) 


with r< = min(rı, r2) and r> = max(rı,r2). As usual, the integral can be broken into two 


segments at rı = r2 


where Ip and Iĝ are 
Ip(v,u) -|/ dry [ drone @ P", (26) 
0 0 
Fe) u’) -f ar f dro rv re our, (27) 
0 


and 


v=a+2+c-q-2k, (28) 
u=b+2+q+2k, (29) 
v =a+2+q+2k, (30) 
u =b+2+c-q-—2k. (31) 
It will be convenient to define the parameter s 
s=v+u+1=v +u +1=a+b+c+5. (32) 


A. Integral Ip(v, u) 


For the case of v > —1 and u > 0, the regular integral Th(v, u) can be expressed in terms 


of the hypergeometric function 


1 OO 
Talos) =i f dry rien RC, u+ 2 Br) 
0 
1 s! B 
- Fı(1 1; 2; ——— 33 
+l lot pps + het 8 a) (33) 


If the conditions v > —1 and u > 0 are not satisfied, J} is singular that can be recast to 


plv, u) z4 eh dra rirse “1 Bro -f im f dro rir3e —arı—pr2 


= Ls(a)L-u(8) — IRV, u). (34) 


When both v and u are negative integers, we should examine the behavior of L_,(a@)L_,(() 


as €a — 0 and eg > 0. A detail analysis is given in Appendix [A] 


B. Integral [},(v,u) 


For the case of v > 0 and u > —1, the integral 7}(v, u) is regular and thus the order of 


the integration can be exchanged, yielding 


s! Q 


mia p | 1;v | eg (35) 


For other cases, see Appendix |B] for details. One can see from Eq. (B5) that Eq. is valid 
for s > 0 and u < 0. Also Eq. is valid for any integer u, provided v > 0 and s > 0. 


I (v, u) = 


8 


For the case of s < 0 and u > 0, the singular integral 1} can be expressed in the form 


I2(v,u) = yee -w+ (+p), @=<Dand a> 0. (36) 


a 
When v < 0, u < 0, and s < 0, we need to consider the singular integrals for some 


special values of (v,u). The following recursion relation can be established, see Eqs. (B7) 


and (B12): 


TRC, =1) = Tail Qi, b) » U= =l (37) 
1 
—u— 1 


where J (x,y) is the auxiliary function introduced by Harris et al. that is given in 


Appendix 


I}(v, u) = [L sla + 8) — BIR(v,u+1)], u< -2 (38) 


IV. DISCUSSION AND CONCLUSION 


Using Mathematica, Iola, b, —2 — n;a,() can be calculated analytically according to 
Eqs. and (18). Examples are given in Eqs. (C4), (C8), and of Appendix [C] for 
Ip(—1, —1, —3; a, 8), Io(0,0,—4;a, 8), and Ip(0,0,—5;a, 8), respectively. Let us consider 
T_50,0(0, a, 8), P_31,1(0, a, 8), and [_411(0, a, 8). T_209(0, a, 8) and '_311(0,a, 8) can be 
directly obtained using Eqs. (13) and (23) of Ref. [IO]. After taking differentiation using 
Eq. (15) of Ref. [10], we obtain 


+ day) 


2 1 2a  3a*- {ork ory py By 
Paleb) = Go { - — + “1h - ng} + ae 
2 L R\2 2 2 _ 2 2 2 
ILS Gy | + OPa + 8+ 812 nla +A) ~ (a+ 9)*Cdory +7? +38") BIn(a+9) 
(39) 
We thus have 
Io(—1, —1, —3; a, B) =T _20,0(0,a, 8), (40) 
To(0,0, —4; a, 8) = T-31,1(0,a, 8) 7 (41) 
To(0,0, —5; a, 8) = T-41,1(0,a, 8) z (42) 


For arbitrary integers a and b, we can use the following formula 


ð a+1 Fa) b+1 
L(G, b, —2 — n; Q, B) — T—1-n,a+1,b+1(0, Q, B) = (-2) (-2) T—1-n,0,0(0, Q, B) , 


(43) 


which can be evaluated analytically by using Mathematica. It is noted that the above 
expression is valid for both regular and singular integrals. 

The formulas for singular integral J,(a,b,c;a,3) with c > —1 have been presented in 
Sec. In order to compare with the formulas for Pim n(a, 8,0), we have derived I9(—2 — 
n, —1,—-l;a, 8), Ipn(—3, —2, —1; a, 8), and Ip(—3, —3, —1; a, 8) in Appendix |C 2] Results of 
Ip(—2 — n, -1, —1; a, 6) with n = 1, 2,3 are listed in Eqs. (C21)-(C23). It is easy to see the 
equivalence between I9(—2 — n, —1, —1; a, 8) and T_j_noo(a, 6,0). It is, however, not so 
obvious to find the equivalence between P_»-19(a, 8,0) and Ip(—3, —2, —1; a, 8), as listed 
in Eq. (C27). Let us first use Eq. (35) of Ref. [LO] 


_ 2 
T_»-10(a, 8,7) =[In(6 + ye — (2)? +1 — ~ dilog (5 7 £) = s 


at+yll,.fat+y ; a+6 
= X Fm (S22) + attog (245) |, (44) 


After taking the limit of y > 0, we obtain 


2 1 
T-2—1,0(a, 8,0) = l+ (1 — yp) + Aa + 5 (Ina + In 8)’ -n’a 


+ 2yg ln 8 — dilog (=£) + dilog (=) + 2[-w(2) +1In 8] lne +ln’e. (45) 


The above equation is the same as Eq. (C27) showing the equivalence of T_2_19(a, 8,0) 
and Ip(—3, —2, —1; a, 8). 


In the derivation of Jo(—3, —3, —1; a, 8), we encounter an unexpected factor €g/€a, see the 
last term of Eq. (C29). This factor comes from the function L_,(a)L_,(G) where v and u are 
both negative, see Appendix [A] The integral [9(—3, —3, —1; a, 8) may have different results 


due to the fact that the expression of €g/€, can be randomly chosen as eg and €a approach 


zero independently. However, To(—3, —3, —1; a, 8) should obey the following relation 


splo(-3,-3, 10,8). (46) 


Moreover, we can determine the expression of €g /€a from the equivalent relation of [9(—3, —3, —1; a, 8) 


and T'_2-20(a, 6,0). Letting y be zero in T_»_20(a, 8,7) of Eq. (B3) of Ref. [LI], we have 


Io(—3, =2, =í; Q, B) S 


T_2-20(a, 8,0) = — (a+ Oh? — 4yr +6 + naln 8 — 4ln(a + B)| — 2(a lna + 81n 8)yer 
+ (8-a) [zma + dilog (=£) | + (a-b) É In? 8 + dilog (=) | 
+ f +2[(2 —yze)(a+ 8) — alna — £ ln £] lne — (a+ 8) ln’ e. (47) 


10 


Comparing Eq. (C29) and Eq. (47), one can see that 


(48) 


which can be applied to other singular integrals involving L_,(a@)L_,(@), see Eq. (34), with 


both v and u negative integers. 


In conclusion, we have developed a general scheme of dealing with the singular integrals 
that appear in variational calculation of three-body Coulomb systems using Hylleraas co- 
ordinates. In this scheme, we first replace the lower limit of an integral by € > 0 and then 
check the behavior of the integral as e + 0. The regular part of the integral does not depend 
on € and the singular part depends on e. We then extend the formulas of Iola, b, —2— n; a, (3) 
from [12] to coincide with Dyimn(a, 6,0) for both the singular and regular parts. Moreover, 
the singular integrals I,(a,b,c;a,() with c > —1 can be evaluated in the framework of 
Yan and Drake [20]. It should be pointed out that the singular integrals I,(a, b,c; a, 3) 
with c > —1 cannot be calculated through Pim nla, 8,0) recursively. We have also ver- 
ified Eq. (6) for Ip(a,b, —2 — n;a, p), Io(—2 — n,-1,-1;a,6), Io(—3,—2, —1; a, 8), and 
Ip(—3, —3, —1; a, 8). In our scheme, all divergent terms that are proportional to either Ine 
or 1/e can simply be dropped, without using a singularity cancellation procedure. It is noted, 
however, there are some singular integrals that should be evaluated using Cy mnn(a, 8,0), such 


as Ip(—3, —1, —2; a, 8) and Ip(—2, —1, —3; a, 8). 
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Appendix A: Function L,(z) 


The general definition of L,(x) is given by Eq. (13). For p < 0, it can be calculated 


according to Eq. (15). For positive integer p, it can be recast to the form 
E, (xe) 
L(x) = a (A1) 
where E,,(z) is the generalized exponential function (see 5.1.12 of Ref. [21]): 


E,(z) = (cay aw = ae i HH , argZ <T. (A2) 


m=0,mén—1 

Thus L(x) can be written as 

(p) —In(we) E_r) __( 00) 
i! Aea 

After taking the limit € + 0*, we can obtain Eq. (14). 
Let us consider Lm(x)Ln(y) which may appear in the evaluation of Ip(v, u), see Eq. (34). 


Lyla) = (-2}°] 


It is easy to obtain an analytical result when at least one of m,n is negative. When both 
m and n are positive integers, some special terms appear in L,,(x)L,,(y) such as €, Ine, and 
€,/éy. In our study, we found that the ratio €,/e, = g(x,y) where g(x,y) is a function of 
x and y. Therefore, special terms of the type €,; ne, = g(x, y)ey, Ine, are deduced to be 
Zero as €s — 0 and e — 0. On the other hand, the expression of €,/¢, remains unclear for 
Lm (x)L£n(y) itself and it can only be determined by properties of referred singular integrals 
I, (a,b,c; a, 8). Subsequently, the explicit formula of €,/¢, is determined as —x/y as €, and 


€y independently approach zero, see Eq. (48). When m = n = 1, we have 


Iy(x)Ly(y) = |Y) — nz -— Ine, — ` = loa —Iny—Ineg, — 5 a | 


k=1 
(A4) 
After taking the limit € + 0 and keeping nonvanishing terms, we obtain 
L(x) L(y) = WY(1) — Ine — Ine,|[(1) — ny — In ey]. (A5) 


When m = 2 and n = 1, we have 


1 = (= Ex 4 = (= Ey j 
La(a)La(y) = (—2)) (2) -Inr — nes — — — 2 loa -iny -lng — 2 me 


= (—x)[y(2) — nz — Ine,|[p(1) — my — ln e] + ty + WA — lny]. (A6) 


For arbitrary m and n, we can evaluate L,,(x)L,(y) in a similar way. 
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Appendix B: Integral [?,(v, u) 


We can perform the integration over rə first for Eq. (27). When u > 0, let us consider 
the integral 


Pe aoe T(u+ 1, 6ri) Ul a, a (Brı)" 
l = [ rye PP drs = Buti = Buti © : os ur Y 20 et) 
$ k=0 i 


where we have used the incomplete gamma function, see 8.4.8 of Ref. [22]. Now I} with 


u > 0 can be expressed as 
ul f° (Br) 
IR(v, u) = Burl f Ti > k! 
G k=0 


l p 
= H1 > k! L_(w+ry(@ +B), u>0. (B2) 
k=0 ` 


e~ (2+8)rı dr, 


For u < 0, the integral [,, is finite and can be expressed in the form 
Ia =r} E (bri), u<0. (B3) 


Inserting [,, into Eq. yields 


I2(v,u) = | rie °"F_.(6ri)dr1, u<O0 (B4) 
0 
with s = v + u + 1. For the case of s > 0, the integral can be reduced to 
Puas aeiee aaa (B5) 
u) = as a) ;>—), s>0Oandu : 
F vt+1(a+ py? a+ B 


When s < 0 and u < 0, the integral I} becomes singular. In practical applications, we 
only need to consider a few values of (v, u). Using the following recursion for E, (see 5.1.14 


of Ref. [21}) 
1 
Enailz) = Ale — zEn(z)], n=1,2,3,... (B6) 


and inserting it into Eq. (B4), we can derive the following formulas for J}. 


In this case, s = —n = —1, —2, —3,... and the integral has been considered in Ref. [10], 


1.€., 


I?(—-n,-1) = (a, 8), n=1,2,3,... (B7) 
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where 


oo et 
Iy(c.y) =f zv(utae, (B8) 


and it can be evaluated recursively according to 


1 2 
I(x, y) = S-F tme- me) +yrllny + Ine) + (Inz + Ine) Iny — dilog (=) : 


2 
(B9) 
pan) = Fe -E lou) - phate +9). (B10) 
2. v=—nu<-2 
When u = —2, I? can be expressed using Eq. 
I2(—n, —2) = ij rêe [e P" — Br, Ey (Bry) dry 
0 
= L_(a+ 8) —BI2(—n,-1), n>0. (B11) 
Applying the recursion relation for E,,(z), we can build a recursion relation for I}, i.e., 
I2(—n,-m—1) = ~[L_,(0 + 8) — BIZ(—n,-—m)], m,n =1,2,3,... (B12) 


Appendix C: Analytical expressions for some special singular integrals 
1. Ip(a, b, —2 — n; a, B) 


In this subsection, we will consider J(a,b, —2 — n; aœ, B) for some special values of a, b, 


and n. For a = b = —1 and n = 1, Oy, O2, and O; in Eq. are 


Ore -20 E, (C1) 
(2 E, "T 
mo eo wo) E 
ani uwe beve 
ERER TE S BO -ijd spe inf-aine), 100 


For a = b = 0 and n = 2, O1, O2, and O; are 


1 (@+8)(ma-Ing) 

O; a (a — 8)? 2(a — B)3 ’ (C5) 
_ v)—mf—Ine pW- inpe. 1 2 

Fat ap (@+B% (ape. (e) 
_YU)-lna-lne afy(2) -lna -— lne] 1 2 

O= Hat) (@+B ttp e. ou 


and thus we have 


4 P AaB i 4aß (a? + B?) 
(a+ Be * FAP A eA) 


Io(0, 0, =4; Q, b) = (In 8 = ii a) : (C8) 


Finally, consider the case of a = b = 0 and n = 3, the corresponding O1, Oz, and O3 are 


B? — a? + 2a p (lna — In 8) 


a 6(a — BP | (C9) 
= p ap oh lie a— Bp | 1 

= asa arp T e age ae O 
= 2 a na+lne— pee 1 

= Sa ae * Her pent me VON t sara Grape O 


and thus we have 


2a 6 


(a? — 62)? (a? f 3a8°) lna — (3° F 3a?°B) In 8 + (a — By lne 
2a8 208 | 2 
~ 3(a—BP(atB) 3(a+ Bee) ECETIA (C12) 


Ip(0, 0, =); Q, p) ~ 3 


2. Ip(—2—1n,b,c;a, 8), Ip(—3, —2,—-1;a, 8), and Ip(—3, —3, —-1; a, 8) 


Let us consider Jo(—2 — n,b,c;a,8) for b = c = —1 and n = 1,2,3. One can find 
that Lə = 0 and subsequently k = 0. Furthermore, the coefficient Ci, is thus reduced to 
C100 = 1. The integral Io of Eq. can be expressed as 


Ip(—2 =n, =]; =]; Q, p) = 2IrR(—2 =n, =l; =l; Q, p; 0, 0) : (C13) 
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Here we list [%(v, u) for some special values of (v, u) according to Eq. (36): 


FR(-1,0) = SU) -nfa + 8) ~ Ind, (C14) 
(2,0) = = - EWO) -mfa + 8) = Ind, (C15) 
[2(—3,0) = sa = — pe ATWO —In(a + 8) — Ine), (C16) 
P(-2,1) = z: : pee Seas) ne T (C17) 
ROB) = sag - ge VO- mate) - ne 

se co E E (C18) 
ROLD = ya a t ope + C WO - me + 5) -mg 

_ (a + BY is Ine hi: (C19) 


Inserting Eq. into Eq. (25), the integral [rz can be written as 
Ip(a, b,c; a, 8; q, k) = L_»(a)L-u(8) — TRl, u) + TRl’, u’). (C20) 


Using the explicit formulas of [?,(v,u) and L,(x) from Eqs. and (14), we obtain the 


following expressions: 


hes 6, B= SWO) nes lama ~(a+B)In(a+f)], C 

fo(4, 1,108) = 542 ~ Ca We) - md} 
+ zle + 8)? In(a + B) — a’ lna], (C22) 

Io(—5, -1,—-1;@, B) = 3E —= | (e + a8 4 =) wa - ng} 

+ ale? Ina — (a + 8)? In(a + 6)). (C23) 

For Ip(—3, —2, —1; a, 8) where s = —1, it can be written as 
Iy(—3, —2, -1; a, 8) = 2[Ip(—2, 0) + I2(-1, -1)] (C24) 

with 

Tp(—2,0) = Lo(a)Lo(8) — Ip(—2, 0) = zlo) = Lala + P)], (C25) 
I(—1,—2) = L(a, 8). (C26) 
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By using Eqs. (A3) and (B9), we finally obtain the following expression 


Ip(—3, —2, —1; a, 8) =1 + (1 — yg) 4 + In 7 > 3 | (na + In 8)? — In? a 
+ 2yg ln 6 — dilog (=£) + dilog (=) + 2[—7)(2) +1n 8] lne + In’. 


(C27) 


Finally we consider To(—3, —3, —1; a, 3) where s = —2. This integral can be expressed in 


terms of L,(«) and I,(2,y): 
Io(—8, —3, —1) = 2[Lo(a)L1(8) — h(a, B) + Lala + 8) — Bhi(a, B)]. (C28) 
By using Eqs. (A3), (AG), (B9), and (B10), Jo(—3, —3, —1; a, 8) can be simplified to 
Iyp(—3, —3, —1; a, 8) = — (a + 8) [ye — 4yr +6 +Inaln@ — 4In(a+ 8)| — 2(a lna + BlnB)yz 
E $ In? a+ dilog (=+) | pñ $ In? 6 + dilog (=) | 


+= +2[2-72)(a +8) -ama - fnf]ne- (a + 6) ne 


+20 +28 8, (C29) 
€ 


Q 


where we have used the following property of the dialogarithm function: 


1 2 
dilog (= r) + dilog (z a r) s= W eL (C30) 
x yY 2 Y 6 
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